function p=predictorMED(a,b,c,d)

%Recibe los pixeles causales a,b,c,d y predice el pixel p de acuerdo al
%modelo MED
%p=median([a,b,a+b-c]);
M=max(a,b);
if c>=M
    p=min(a,b);
elseif c<=min(a,b)
    p=M;
else
    p=a+b-c;
end